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RESEARCH ACCOMPLISHMENT IN THE FIRST YEAR 


The research program was planned in four stages over two calendar years 
starting September 28, 1979. The first and second stages would encompass both 
the investigation of free surface properties of the semiconductors in an altra- 
high vacuum (lO"^” torr) system (UHV) equipped with LEED, Auger electron spectro¬ 
scopy (AES) and ion-sputtering capabilties, and the design and assembly of the 
satellite chamber for device fabrication. These two stages would be performed 
at the same time without interference with each other. The third stage of the 
work would encompass the study of the waveguide surface structure, formulation of 
surface-scattering theory and characterization of scattering loss parameters in 
terras of the surface geometry through LEED and SEM analysis in the same ultrahigh 
vacuum system. The last stage of the work would encompass the fabrication of the 
oxide-free electro-optical devices, the measuring of the optical propagation 
losses and the study of the optical characteristics of such waveguides. 

The first and the second stages of the planned research were completed. 
Currently we are conducting the third stage of the research plan. The satellite 
chamber for device fabrication, as shown in Fig. 1, was designed and assembled 
to the main system. A schematic diagram of the main system is shown in Fig. 2. 

The satellite chamber is indicated as subsystem 2 in Fig. 2. Subsystem 1 in Fig. 2 
is a plasma reactor designed for oxidation research of GaAs and is in no inter¬ 
ference with this research program. A photograph of the completed system is shown 
in Fig. 3. 

The free surface properties of GaAs (110), GaP (110) have been studied by LEED 
and AES. Results of the studies have been presented in two different scientific 
conferences. Abstracts, which summarized the results of the research works, are 
given in Appendices A and B. Abstract of the first paper entitled "Chain Method 
of LEED/MEED Intensity Calculation for Diatomic Surface" is given in Appendix A 
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SPECIMEN PLUG 


Fig. 1. Top-view schematic of the satellite system attached to the 
main chamber (left). A: sorption forepump; B: Vacion UHV pump; 

C; bellows-coupled flange; D: Isolation valve; E: six-armed cross 
F: ionization gauge; G: transfer probe magnet actuator. In the 
schematic the transfer probe is shown completely retracted, with the 
specimen plug located at the center of the six-armed cross. 
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1^1 g. 2. Schematic diagram of the UHV system equipped with LEED, AES, mass 

spectrometer, ion sputtering gun and ellipsometer for surface analysis 
Subsystem 1 is the plasma reactor which is currently under constructio 
Subsystem 2 is the ‘CS device fabrication chamber. constructio 
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LEED/AES system for Investigation of surface ami interface properties of semiconductors 
and ceramic materials. Tl>e capital equipment is acquired througi» Uie USE equipment 














































and the second paper entitled "Atomic Structure of GaAs (110) Face" is given in 
Appendix B. Both papers are being prepared for publication. 


The third stage of the research work is currently ongoing. The theoretical 
formulation of surface scattering theory and characterization of the scattering 
loss parameters are currently being developed. The problem of a wave propagating 
along z direction in a rectangular dielectric waveguide with dimensions a and b as 
shown in Fig. 4 was formulated and solved. The wave functions in the waveguide 
were solved using predominantly polarized approximation.* The propagation 
attenuation constants, e.g. the loss parameters, were calculated. Details of the 
theoretical formulation and calculation is given in Appendix C. 

Citations of research results, which are presented in conferences or sub¬ 
mitted for publication, are listed in Appendix D. 


i 


*D. Marcuse, Bell Syst. Tech. J., 3187 (1969). 






Fig. Cross section of a rectangular dielectric waveguide 
with dimensions a and b. n], ri2> n3. n4 and nc are 
the indixes of refraction of regions 1 to 5, respectively. 
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APPENDIX A 


The following paper was presented in the Conference on Determination of 
Surface Structure by LEED, ISM Thomas J. Watson Research Center, Yorktown Heights, 
N.Y., June 19-20, 1980 

CHAIN METHOD OF LEED/MEED INTENSITY CALCULATION FOR DIATOMIC SURFACES 

In low energy electron diffraction (LEED), crystal surface is probed by bom¬ 
barding electrons of energy E £ 200 eV at normal incidence to the surface and analy¬ 
zing the intensity of those elastically scattered electrons which are back reflected 
from the surface, whereas in medium energy electron diffraction (MEED) due to an in¬ 
crease in the energy of incident electrons (200 eV < E £ 5 keV) the surface sensiti¬ 
vity of the technique is maintained by taking an oblique angle of incidence. Inter¬ 
pretation of LEED/MEED intensity spectra requires detailed calculations of the 
diffracted intensities for a series of trial models. For sucn a technique to work, 
a rapia and accurate method of calculation is needed. 

Layer-KKR method is conventionally employed in the calculation of LEED intensity 
spectra. In this formalism, the crystal is divided into a number of layers parallel 
to the crystal surface and scattering calculations are split into two parts, intra¬ 
layer and interlayer. In the intralayer scattering calculation, due to the assumed 
spherical symmetry of the atomic potential, an angular momentum representation is 
used. This involves large matrices and their inverses and becomes cumbersome with 
increasing energy and more complex surfaces. In such a situation we propose to use 
the chain method of intralayer multiple scattering calculation. 

In the chain method the two dimensional intralayer multiple scattering calcu¬ 
lations are further divided into two one dimensional steps, scattering within a 
chain of atoms and between the chains. At each stage of calculation the electron 
wave function is represented in terms of an appropriate set of basis functions. The 
scattering by an atom, a chain of atoms and a layer of chains is expressed in terms 
of spherical, cylindrical and plane wave representation respectively. As we move from 
one stage to another, transformation from one basis set to another is carried out. 

Due to its one dimensional lattice summations the chain scattering formalism has 
many computational advantages over the layer-KKR method both for normial incidence 
LEED and off normal incidence MEED. 




APPENDIX B 


The following paper is submitted for presentation in the Eighth Annual Con¬ 
ference on the Physics of Compound Semiconductor Interfaces, Williamsburg, Virginia, 
January 27-29, 1981. 


ATOMIC STRUCTURE OF GaP (110) FACE 
ABSTRACT 

Low energy electron diffraction (LEED) intensities have been measured for the 
(110) face of GaP and analyzed using a dynamical multiple-scattering model of the 
diffraction process. The intralayer multiple scattering is treated exactly, while 
for the interlayer multiple scattering, the renormalized-forward-scattering method 
is used. Comparison of the calculated and observed LEED intensities suggests that 
both the Ga and the P atoms on the (110) face may exhibit a contracted outermost 
layer spacing. The surface l^y^C ''s compressed by about 5% such that the top layer 
spacing is reduced by 0.1 ± 0.02A. The rippled geometry of surface reconstruction 
is not clearly observed. This indicates that the GaP (110) surface atomic struc¬ 
ture is different to that of GaAs (110). By comparing the structures and properties 
of the GaP (110) face and the GaAs (110) face, it is concluded that the GaP (110) 
face is relatively unstable and reactive. 









APPENDIX C 


The problem of a wave propagating along z-direction is a rectangular dielectric 
waveguide with dimension a and b as shown in Figure 4 is formulated and solved. An 
exact analytical treatment of rectangular waveguide is practically impossible. 
Therefore the approximated analytical approach developed by Marcatili (a) is followed 
in solving the problem. 

Assuming a traveling wave propagating in z-direction as shown in Fig. 4, 

Maxwell equations 


7XH = II 




can be solved to give 
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where n is the index of refraction of the dielectric; k is the wavevector in 
free space; oi is the angular frequency of the wave and 3 is the propagation 
constant in z-direction. 

Substituting (2) and (4) into (5), and (1) and (3) into (6), we have 


3^E 3^E 

=0 .(7) 
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Two different types of modes can be supported in the waveguide^ , e.g., 


Epq : polarized predominately in x-direction and 



: polarized predominately in y-direction, 

HH 

where p and q are positive integers indicating the modes of the wave propagating 
in the waveguide. 

eL mode 
pq 

In rogion 1 with refraction index n^ as shown in Fig. 4, the waveequation (7) 
can be solved to give 


E 2 (x,y) = A Cos {x+ 5 ) Cos (y+n) 


with = 0 , where C and n are the phase factors of E^. 

Other components of the wavefunction can be obtained by substituting {9a) 
into eqs. (1), (3) and (4): 


H = - Sin k„ (x + E) Sin k^, (c + n) 


*^o X 


(n?k^-k^ ) 

Ex = ■'A - — Sin k^ (x + E) Cos k^ (y + n) 


Ey = —Cos k^ (x + 5 ) Sin k^ (y + n) 


Substituting (9a) into (7) we have 
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For small incident angle, S»(k^ + kM, eq. (10) can be approximated to be 

X y 

3 = n^.k Cos0 = n^k for i<S°. 

The wavefunctions in regions 2, 3, 4 and 5 with indices of refraction n^, 
n^, n^ and n^, respectively, as shown in Fig. 4, are solved. When boundary con¬ 
ditions are matched between region 1 and the neighboring regions, i.e., regions 2, 
3, 4 and 5, we obtain the four transdental equations; 


tan k^a = n^k^ (n^yg + n^Y 3 )/(n|n| kj - n:[Y 3 Y 5 ) 

tan k| = (k^/Yj). 

tan kyb = ky (y 2 + Y 4 )/(k‘ - Y 2 Y 4 ) . 

tan k^n = -Y^/k^ . 

where Y 2 = C(n| - n^)k^ - 

V3 = [(nf - n|)k^ - 


(lla) 

(lib) 

(12a) 

(12b) 


Y5 = [(n? - n5^)k^ - 

A computer program is currently being developed using Newton's method to solve 

equations (11) and (12) to obtain k , k , 6, and the power attenuation constant a. 

X y 

REFERENCES: 

(a) Marcatili, E. A. J., Bell Syst., Tech. J. 48, 2071-2102 (1969). 

(b) 0. Marcuse, "Theory of Dielectric Optical Waveguide", Academic Press, New 
York, 1974. 

(c) I. P. Kaminow, W. L. Mammel, and H. P. Weber, Applied Optics, Vol 13, No. 2/ 
Feb. 1974. 
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APPENDIX D 


Publication citations and presentations supported or partially supported by 
this research grant: 

1. "Decomposition of Aluminum oxide by Electron Bombardment." B. W. Lee and 
J. M. Kuo, B. Am. Phys. S ., 25 (3), 238, 1980. 

2. "Study of MIS Polycrystalline Silicon Solar Cell Using Auger Electron 
Spectroscopy," J. M. Kuo, B. W. Lee, B. Lalevic and W. A. Anderson, B. Am. Phys. S ., 
^ (3), 409, 1980. 

3. "Chain Method of LEED/MEED Intensity Calculation for Diatomic Surfaces," 

N. Masud, C. G. Kinniburgh, 0. J. Titterington. Presented in Conference on Deter¬ 
mination of Surface Structure by LEEO, IBM, T. J. Watson Research Center, 

Yorktown Mights, NY, June, 1980. 

4. "Study of MIS Silicon Cell by ESCA and AES," Y. S. Wang, H. J. Yu, C. C. Hsu 
B. W. Lee and W. A. Anderson, presented in the 27th National Symposium of American 
Vacuum Society, Detroit, October, 1980. Submitted to J. Vac. Sci. Technol for 
publication. 

5. "Stability of MIS Silicon Solar Cell," B. W. Lee, J. M. Kuo, B. Lalevic and 
W. A. Anderson, submitted to J. Vac. Sci. Technol . for publication. 

6. "Atomic Structure of GaP(llO) Face," B. W. Lee, R. K. Ni and N. Masud, 
submitted for presentation in the Eighth Annual Conference on the Physics of 
Compound Semiconductor Interfaces, Williamsburg, Virginia, Jan. 1981. 
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A NONLINEAR MAXIMUM ENTROPY METHOD FOR SPECTRAL ESTirvATION 

Summary 

An intensive research study is made of the nonlinear maximum 
entropy spectral analysis method proposed by P. F. Fougere of 
AFGL. The research not only provides a much better understanaing 
of the properties of the method including spectral resolution, 
convergence, etc. but also shows that the method can be perform'.'! 
effectively with the PDP 11/45 minicomputer. The successful 
implementation of tne method is described in the program listings 
in Appendix I. Extensive computer results are presented on 
various data. These results clearly confirm that the nonlinear 
method is superior to the Burg's maximum entropy spectral analysis. 

Preliminary results on the multichannel (multivariate) maximum, 
entropy spectral analysis, the two-dimensional maximum entropy 
spectral analysis, and computer graphics for the spectral display 
are also presented. An extensive bibliography of the miaximum, 
entropy spectral analysis is given in Appendix III. 
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I. Review of Research Progress 

This research is concerned with the nonlinear maximu;n entropy- 
spectral analysis (MESA) method proposed by Dr. Paul Fougere of 
AFGL. As verified in this research work performed at the 
PDF 11/45 minicomputer, the method not only provides a much better 
spectral resolution than the Burg’s method but also removes the 
line-splitting and frequency shifting phenomena for sinusoidal 
signals, as experienced in the Burg’s method. By using tr.e double 
precision, the minicomputer results are reasonably close to those 
obtained at the GDC 6600 computer at AFGL. 

The final computer program developed for the nonlinear 
complex sigrial maximum entropy spectral analysis is chowr. in 
Appendix lA. The program follows the mathematical development of 
Fougere [ 1 ] but is quite different from the original computer pro¬ 
gram developed by fr. Fougere. Appendix IB is the computer 
program for Burg's complex signal maximum entropy spectral analysis. 
By using the computer programs and the two-channel radar aata as 
shown in Fig. 1 (see also[”2] ), the spectrum of the nonlinear 
method is shown in Fig. 2a (linear plot) and Fig. 2b (logarithmic 
plot) for 10 filter weights. The Burg's result is shown in Fig. 3a- 
(linear plot) and Fig. Jb (logarithmic plot), also for 10 filter 
weights. More detailed tabulation of the major frequency compo¬ 
nents outside the clutter bandwidth (-0.l67fs, +0.l67fs) is given 
in Fig. 4. The nonlinear method which is very close to the true 
answer clearly is much better than the Burg's method. 30 iterations 
are used in the nonlinear method which appears to be an optimum 
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number. The optimum filter weight is around 10 or 11 as 
determined by Fig. 5 which shows the linear prediction error (FiZ) 
according to the Akaike criterion. The solid curve is the lower 
bound and the dashed curve is the upper bound as tabulated in 
Fig. 6. The upper bound was originally proposed in the research 
proposal [2] and it seems to be better than the lower hound due to 
Akaike. 

Major documentations already made which describe the research 
progress are as follows! 

1. C.H. Chen, J. Chen, and C. Yen, "A minicomputer implementation 
of Fougere's maximum entropy spectral analysis method," Techni¬ 
cal report prepared for the Mini-grant, August 20 1980. This 
report has detailed results of comparison between Burg's and 
the nonlinear methods for sinewave and sunspot data. 

2. C. H. Chen, "Spectral resolution of Fougere's maximum entropy 
spectral analysis," to be published in the Proceedings of IEEE, 
June 1981. This journal paper based on the work performed 
under the Mini-grant provides a good comparison between the 
nonlinear and Burg's methods, and the other method for complex 
sinusoids. The Cramer-Rao bound is used as a reference. 

3. C. H. Chen and C. Yen, "Note on computer graphics for maximum 
entropy spectral analysis," Technical report prepared for the 
Mini-grant, March 23» 1981. This report provides a three- 
dimensional spectral display of sinewaves for both nonlinear 
and Burg's methods. 
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A number of important results are included in the new research 
proposal submitted to AFOSR in December 1980. Appendix III 
provides an extensive list of references on the maximum entropy 
spectral analysis. The following sections describe some new 
research areas with preliminary results. 

II. Multichannel (multivariate) Maximum Entropy Spectral Analysis 

A mathematical presentation of this topic is given in Appendix 

II. Several computer programs for multichannel maximum entropy 
spectral analysis were provided by Dr. Fougere. The following 
results are based on the time series of sunspot numbers, northern 
light activity, and earthquake activity by using the third multi¬ 
channel program. The data are tabulated in [3]- 

Fig. 7a is the first channel (sunspot number) auto-spectrum 
with linear (left) and logarithmic (right) scales, and l6 lags. 

Fig. 7b is the second channel (northern light activity) aumo- 
spectrum with linear(left) and logarithmic (right) scales, and 
16 lags. 

Fig. 7c is the third channel (eartnquake activity) auto¬ 
spectrum v/ith linear (left) and logarithmic (right) scales and 
l6 lags. 

Fig. 7d is the cross-spectrum between channels 1 and 1 with 
real part (left) and imaginary part (right), and l6 lags. 

Fig. 7e is the cross-spectrum between channels 2 arid 3 v/ith 
real part (left) and imaginary part (righ ), and 16 lags. 

Fig. 71 is the cross-spectrum between channels 1 and 3 v/ith 
real part (left) and imaginary part (right), and l6 lags. 

Dy way of verification, it is interesting to note that the 
spectral peaK for the sunspot numbers is determined accurately. 

III. Two-Dimensional Maximum Entropy Spectral Analysis 

he consider a very simple separable case in the tv/s-dimensional 
spectial analysis. The signal considered is sin(2Fx) sin(2?ry). 

In tnis case the power spectrum is the product of the power spectra 
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of sin(2T[x) and sin(2Try). In each spatial dimension, the power 
spectrum can be determined by using the Fourier, nonlinear and 
Burg's methods. The two-dimensional spectra are shown in 
Figs. 8, 9 and 10 respectively, based on the Fourier, nonlinear, 
and Burg's methods. The nonlinear method clearly is much better. 
Extension of the above procedure to a more general two-dimensional 
spectral analysis is not possible. Although some two-dimensional 
maximum entropy spectral analysis work has been reported (see 
Appendix III), the success is very limited. Further research is 
much needed. 

IV. Conclusions and Recommendations 

The nonlinear maximum entropy spectral estimation method 
proposed by P. F. Fougere has provided superior spectral estima¬ 
tion over the Burg's method in a number of data considered. The 
computational requirement of the nonlinear method is, however, 
significantly higher. Typical number of iterations needed is 20 
to 30 . The use of minicomputer has not created significant 
computational problem as predicted. It is our strong belief that 
the nonlinear method will become very popular in high resolution 
spectral analysis for a wide range of applications in geophysics, 
sonar, radar areas, etc. 

Recommended further studies include the multichannel 
(multivariate) maximum entropy spectral analysis, the two- 
dimensional nonlinear maximum entropy spectral analysis, computer 
graphics for the spectral display, signal decomposition, and 






signal prediction and extrapolation. 
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No. of weight 

Log. of 

n+Tm+i) _ 

Log. of 

N+2(M+1) 

N-(M+1) 

N-2(M+1) 

1 

-0.279 

-0.224 

2 

-0.463 

-0.381 

3 

-0.445 

-0.332 

4 

-0.562 

-0.418 

5 

-0.674 

-0.496 

6 

-1.080 

-0.866 

7 

-1.539 

-1.283 

8 

-2.063 

-1.762 

9 

-2.587 

-2.006 

10 

-2.597 

-2.175 

11 

-2.678 

-2.176 

12 

-2.686 

-2.075 

13 

-2.738 

-1.971 

14 

-2.815 

-1.764 

15 

-2.839 

- 

16 

-2.824 

- 

17 

-2.833 

- 

18 

- 2.810 

- 

19 

-2.799 

■ 


Lower Bound Upper Bound 


Note: N is the number of data points in each channel (32 in this case) 
M is the number of filter weights 
is the error power 
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i’^pondix IE Bur-j's Coiiplex Signal. tlRSA Conputer Program 


■OKI KAN 

00’>:.. 

1.3 (.1(3 :.-;(.i;j 4 l'5-l'iAK-8 ,l rAi.4- 


C: 

NOI ih : ChPBIJiA C :>. 1, MI DA i> ; 1 'I'-AI .IG-80 

>001 


iOl'! IN)Hr X 

■ *oo:? 


nKAI. P ( 3:;') , Pbh (32), PFK ( 32 ), FNk W ( 32 >, GMh 0 < 32 ) 

1003 


rvhM-'-il U () j U'j'L:’ ( ) 1 1 j*Ii ( ) < HAl.f"' ( ) ■ H (i i'.' ) 

•004 


KhA!, Sh (513), SOflS:., YX( IdOO), CA!.( : OOo) 

1 >00b 


UU! li-'l h.X hi Pt: h j PLKj LiLil.P l.iU'.iil» P; Hf’il ..h'/ H 

•003 


UK I 1 h ( 6, O') 

;)007 


FL-KNAT ( 1 X, ■■•.(■■>)■It- I.'I-' HlJli. LG. 'rNj bXPAMii, )J > 

>003 


KPAiKA, lOMIUn, 1 G, '31.!, PXPAIIU, DT 

'.<009 

10 

hi iPnA l (215, 3P;.?. 5) 

■:>o 0 


WKJ I r. (6, 15) 

.'01 ;i. 

15 

rOKMAl (13, -if-I/-' I'dIAKI, rOiUP, I h'l., IIU J 1 , 12 

(1012 


RF AD ( 6, 20 ) 131 AP I , 13 rop, NOl 11 , 5?, I.P 

>013 

20 

Fl'Pi'ir'i l' (6)5) 

.'Ol'l 


NK.,: iP'b, 25)NUM, LG, IRTAK ), P-. TOP, INC, NGU : , I.A 'ON 

>015 

J. 1' 

P.:.KliA rMHl, 5X, /17, OR 12 . 

■-<013 


CALL. BFI.l 

(>01 7 


N' r 'jl.iLKP > 

■1013 


DPR I nr: R ,1 LE 2 (, MI, LI IN)»- X ) 

0(319 


CALL. INPUT ( R , NUM 3N. RNRL ■ 1 .-NR.II) 

.1020 


WKT rR.(5, 3(3)NIJM 

0021 

30 

i" i.tpi'iA f * F EP'i 1 A 1 i.*!”' h’1“ ; 1‘rv) 

0022 


CAI.L. WRlNUn, R, 5) 

0023 


I PI'^8. -(A l AN( 1. ) 

'3(324 


UAi 1... RiRI.,.! 

0025 

5(3 

CON r ;i NOP 

0026 


D R! •.) ■■■=> R1 - 0 A ■ 1 ( M C U 1 - 1 ) * R. X P A1' 1 > 

0027 


RRlJlu dN 1 (LiRN./H;XPAND--1. ■ 

(3028 


D-' 60 Kr >' 1, PRUI 

( . (329 


3R (')•• R ) ~R D.iA 1 ( PI— 1 )./I.iR N—. 5 

0030 

60 

3R (R'RO-rp] (^'SRdPR ) 

0031 


DO 65 I>-/,NLIM 

' (J032 


RPRIH I )-=KR Al... (R ( I ) ) 

• .3033 


CNRIJ ( I •"•=A THAI PRO 1 ) ' 

0034 

6 5 

C( >M r ] N JR. 

! 0035 


DO 70 I-'>,NOri 

(3 (3 3''6 

?0 

YX (I I'^RLOAT ( I :■ 

I 003/ 


CAI 1. NRNPA''l 

0038 


CAl L LDPI Ul (YX, PNRl'J, RIUN, 50, lOKI, 5(3, 'TOO' 

(3(;K-'9 


CAI L R IGl ( YX, CAL, 3-1E. 3) 

(3040 


CAI 1 BhI 1 

i 00^1 


Rip I 1 E ( 6, ) 

(30 <(12 


ROKRAT ( 2.:iv, I./p CHi'NMRI. i NAVE -oh'*-) 

0043 


KEAl Gc', '->'-'9)NX 

00^ 4 

vv V 

Ro 'sMA 1 ( 15 :■ 

CiO^I 5 


oALI. NKUPAo 

(3(346 


CALL PPPI I 'l (YX, CNRI-- NUi'L Ro, lOiO, LO:. , /CO) 

(304 / 


CAL L R IGl ( YX, ( AL, ::43, 3) 

0048 


CALL. DPI 1 

0049 


NKI I F (6, ■■> ■0 ) 

>3050 

9 7/' 

RGkNA 1(20/, 1 /-■ CHANNEI T[ WAL'R: iOHH/(■ ■' 

0(351 


READ (6, ■^■-'■-')LL 

(3(352 


K -0. 

(3(3 5 .'5 


JJU R>0 I - 1 , MOP 

i 0(354 

80 

p”,-'-i R ( 1 ) itiioN.jioa- ( I )) 


/ 
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ORTRAN 

VOo. J:-: 

00 30; 14 25-MAf<-8- 



P'^'P/ELCA'l (NI.IM) 



8AVE'=P 

0':>7 


WRITE (5, ;E5)p 

'058 

•li*!!' 

1- oklia r (21X, •• p - , r 1 5 7 ) 

'.<59 


LAI-MOM-1 

" lAC) 


DO 150 MIL-i,LA;i 

■001 


nnm-nm+i 

'062 


CALL.. Oh‘KL'(h!UMi PhL# PhHi GGUi’^ij H) 

'0/-. 


P- < 1. -G (IvllM ) -HOuiM.JG (>:- ( Miv ) ) ) ■'rp 

'064 


HLI.P-P 

065 

:i 00 

IFiMI'IN Cl .18 it Ip) CO TO 150 

v-66 


IPdVNM. L'l i 81 AP I )CO TO 15 0 

'067 


lEiMGLK (NiJN-IS'lAR'iILIC). ME. 0)00 Oj ; SO 

06:-: 


Lfil-L oPh( Nl'j 1 lij Li oH i KLUl 1 P- D \ ) 

069 


ChI 1. EEI...I 

070 


Col 1... ):-;l-l...l 

■071 


RL Ai) ( 6, 'r'9 7 ;■ MX 

'072 


CAI 1. MEWRAC 

■0/3 


CALI. KLiPl.i I r ( YX, 8, NnIJT, 50, 1010, 50, '00) 

10/4 


WRI TEO-., l::'0)MMM 

■■;)V5 


WRllECj, 120)MtV 

'076 

120 

hL*f‘<r''Al (L^OXi *'-i;--i'r«-if i\ 11 Qp [3. 

'077 


CAI L MR (MM, C, 15 ' 

•078 


CAi L E IGl ( YX, CAL, 1 2, 2) 

»079 


IE (IP. EO. 1, ) CAI. 1., PR IT (MOU i ■ D't, S, Y X :• 

^Oc!() 


IE (EC. LE. 0)G0 TO 140 

1081 


Rf-AD<6, 9'-'9)MX 

'0:52 


CALI ALC(8, MQl.l 1 ) 

.'083 


CALL MhMRAi-i 

.>0:84 


CAIJ. WR(NM,G,5) 

085 


CAi 1. KBPL.01 < YX, 8, MiiUT, 1:0, 1010, 50, ^OO > 

>0:86 


MRITE (6, 120)I'.!MM 

>08 7 


CALI FIGl (YX, CAE, 12, 1 ) 

.'OS 8 

140 

E-HELP 

>089 

150 

i Oi'l 1 INI IE 

>090 


Col L I3ELI 

>091 


CAI..L EXI i 

>(>92 


Ei'E.i 


Kill IT INK 

8 CALI.Elj; 


HI-I..L , 

IMEUT , MR , AT AM , El J lAi , .iMi 


AI MAG , 

IMEMRAG, LBRLOT, EIGl , Ci'NJG , BPEC 


SPEC , 

PPM , AI G , EXIT 


C'P T :i 0N:5 

= ./ijP: 2 


BLOCK 

L.EMG 1 .H 


MAIM. 

8185 (03V762)<^ 


. $4$$. 

2 (000004) 


■«--»CCi|-'|HlLEh'-COKI-^H'- 

PHAsh iJstD l-Pl-E 

DPI :LAKAr JVrS OOS-SA 14'51? 
EXEi :n rAn:|..KS 01263 1411':- 
ASPhliBLY 01 '74 18244 
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t . 


1 

j 

i 


i 

i 

1 

■< 
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' ' 'F<TKAN VGA. 


' >01 
' )0'/ 

' >03 
>04 
■ )0b 
>06 
> 0 / 

>03 

■i.iOO 

010 

till b 

>01 2 
■.-13 

■Ml 4 10 

•lib 

■016 20 

■01 / 

■Ml 

'01’? 30 

■■.).<-0 40 

021 
‘022 


CO ::i 4 4 2b-l'AK-3l 

SU).-:Ki:iU l INK INHUT (!<, HI , SH. Y, L ) 

CiOl'-VOOlM INIltX 
COl'ir'LhX XCi) 

KtOL V(N■),/(N-) 

INOhXtt:! 

RK0;:'(2 ■■■ IIVDF-X )Y 
IlviUhXtt:2 

6600(2 ■■ I lv‘))h.X ) Z 
DO 0=1, Ml 
Z6!- . —Z( J) 

COM r INI 16 
DO 10 0 ■ 1 ■ F' 1 

X( I )=CM6LX ( Y ( I ), Z ( I ) )+SIM-;h:M6I X (GAUCSC 
CUM 1 II\1IJ6 
Wh' 1 i fc ( Ij I ..''..I) 

FOKi'lAT (/'.'MX, CHAIMNi-L I !-X. ' CHAIMVl-i T 
DC.! 40 ,i -1 , M i 
U 6 n £ (, IC )/(!), Z ( I ) 

FOKi'iAl ( " C X 6V, ;i, C'X, K’:-'. 6i) 

CCMIIMOh 

RHI'IJRIV 

FMU 


ROI.11 ;i M£S i.;:a; .l.kd: 

CM6LX , GAI.'S3 

OPT III NS ■- / OP: 2 

Bl OCl, L.6NG I 61 

IMPlJI 2/:: (0OiO42>->t- 

. 2 (OOi'OO^ ) 


»^tCONPIL£P — 
PHAS6 

DECI APATlVhS 
b XbCO I AfiLhS 
ASSEMBI .Y 


- C0R6^ii--:'r 

036n 6666 
0062.6' 14 Zb'.' 
0('7S6 14f.M2 

012 IS 13S00 


1 

i 



/ 


6A> ’6 


1 . i (.■ >, OhL 

I • ■ 
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■ RTKAN 

VOX. lb 

00; 32; 11 2':>--l'iAf- 

1)01 


;3I.IEKi;iU ri,Mr NK^NNr O, IC) 

.)02 


Rb AL O( 1 ) 

OOb! 


CGNPI EX i3 

> 04 


DO 1 0 I E NN 

OO'b 


A-REAI (G< I ) ) 

>06 


)-:-AII'iAi:'i(0( f ) ' 

007 


[b ( IC. EO. '.5)WRnb(5. 20)A, E 

OOP 


IbdC. EQ, A)UKr!E(o. 20)A, B 

009 

20 

bOKNAT ( 2:3X, ■" ( ", b 15. 7, ', ■', b i: v 7 . ■ ) 

0 3 () 

10 

CONIINUb 

oil 


RETORN 

03 


END 


ROUlINE 

S CAU.ElJ; 


KEC^L , 

AI NAG 


Ob'i ION:-: 

^-./OP; 2 


El ( II :K 

LENG'I H 


WR 

3 50 (000454 )'*■ 


irttCOhl-'ILEK - r, f IRh 

PHASh U3KD hRhE 

HKCLARAT) V-S 0070 '.* 
b Xh CU I AOI HS 00 XP3 1 
ASSbMhU.Y 0109 V 13921 


bORTRAN 

V06. 3 3 

(>0; 32; 31 

'■>00 :i 


bUNCTICN GAU::;S( XBAR. 3IGI3A) 

1)002 


R00 1 = l ./Pi'RlOr yATANCl. )) 

ooo;7' 

10 

corn' II lOb 

1)004 


< -30 ■)■• ( Rl'':i'U ( >1 1 — 5) 

0005 


Y -'’;i lOT-'-fiXP (->dX./2 :■ 

0006 


YI RY-4vi:ii:'i'l -tt-RANd). 3 ) 

0007 


Ib(YIRY iVI Y)O0 lO 10 

ooos 


OA; i:-.9=X#:t:IGNA-i-X)::AR 

0009 


Rb TURN 

OOlo 


END 


ROUfJIMbS LAIJ.hD: 

SORT , ATOM .. RAN , bXP 


iipriLNP. -'OP 2 


El HOP l,t:;NGiH 

iiAUSP I'l \ oOn^: Z-R) * 
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Appendix II 

On Multichannel (Multivariate) Maximum Entropy Spectral Analysis 


1. Introduction 

The univariate maximum entropy spectral analysis has now been well developed 
and applied to many defense research areas such as radar, sonar and geophysics. There 
has been some work done to extend the maximum entropy spectral analysis to multivariate 
case. Whittle [1] and Robinson [2][3] generalized the Levinson-Durbin recursion to 
the multivariate case by fitting both foward and backward autoregressions in a 
stepwise fashion. In this thesis. Burg [A] has mentioned about the multichannel case. 
However the computer programs for both multichannel and multivariate maximum entropy 
spectral analysis were only recently developed successfully. Morf et.al [5] developed 
an algorithm for direct estimation of the normalized reflection coefficients from the 
observed data for maximum entropy spectral analysis. They also compared the spectral 
estimation with the methods of Jones [6], Nuttall [7] and Strand [8], which are more 
of a direct extension of Burg's work to the multichannel (multivariate) case. Burg's 
algorithm does not generalize directly since the forward and backward autoregression 
matrices are not the same in the multivariate case, and the forward and backward one- 
step prediction error covariance matrices are different, although they have the same 
determinant. In this report, the programs developed by Strand and Jones are applied 
to real multichannel data and imagery data in addition to a set of test signals. 

The merits of these methods are closely examined. In spite of programming complexity 
the multichannel and multivariate maximum entropy spectral analysis will have increased 
application as the real data are almost always gathered in several channels. Data 
from several channels form a vector for multivariate study. 

2. Brief Mathematical Analysis 

Let x^,X2 ,...,x^ denote n zero mean vectors of dimension d each. The sample 
estimate of covariance sequence for lag j is 


/ 
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1 

2 n t+j t 


( 1 ) 


where the prime denotes the transposed vector. The forward and backward predicting 
autoregressions of order p are given, respectively, as 

^t-k 


= I A 
k=l 


-Cb) = 


' ' = I \ 

^ k»l ^ 


( 2 ) 


t+k 


where and are d x d matrices, and can be determined recursively [6] by 


making use of the estimated covariance matrix in Eq. (1), The recursion starts with 

o o o 

The one-step forward and backward prediction error covariance matrices are 


(3) 


= (I - 

p P P p-1 

s'") . (I - B'P)A'Pts«>) 

P P P p-1 

The forward and backward residuals are, respectively, 


(A) 


x^ - I A^P^ X 
t ,k t-k 
k=l 


t = p+1,...,n 


(5) 


^t^^ \ J ^k^^^t+k’ t-l,...,np 

k=l 

The recursive equations are then given by 


.(p) 
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n(p-l) 

®8-p ’ 

t = p+1,...,n 
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least 

squares 

estimates for the 

forward and backward 

autoregression matrices 
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= uv-^ 
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3(P) 

= u'w“^ 






where U is the sum of cross products of forward and backward residuals at lag p, 

tf'"-”)' (8) 


t=i 
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and V and W are estimates of (n - , (n - respectively, 


n-p 


V = I B, 

t=l 


(p-1) 

t ' 


(9) 


" - T 

t+p t+p 


( 10 ) 


Although the forward and backward autoregression matrices and the prediction error 
covariance matrices are different, the multivariate spectra should be identical 
when calculated from the forward and backward fits by 


S(f) = h[A(f)]“^ [A*(f)] 

or by 

S(£) = h[B(f) ]"^ [B*(f) ] 

where 

A(f) = I - I A^P) 
k=l 


B(f) = I 



-2TTikhf 

e 


h is the sampling period and * denotes complex conjugate transpose. 

The above approach based on the work of Jones [6] does not guarantee stability 
and does not generally produce a non-negative definite spectrum as has been pointed 
out by Nuttall [7]. Subsequently, Nuttall [7] and Strand [8] applied a weighted 
arithmetic mean error criterion in order to provide model stability and to ensure 
positive definite stationary spectra. Another procedure suggested by Morf, et.al. 
[5] that also meets the spectral requirements is to compute the spectrum from the 
normalized reflection coefficient matrix p. To obtain this matrix, W and V are 
factored by using Cholesky decomposition into the product of lower triangular 
matrices times their transposes. A new recursive procedure for and by 

using p in place of Eqs. (A) can then be obtained. Other recursive algorithm has 
been proposed [9] for the solution of the normal equations for both single and 



/ 



multichannel data. 
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Appe’^.dix III 


A BIBLIOGRAPHY ON I.iAXIKUM LNTRCPY BPLCTHAL 
ANALYSIS AND HLLATSD TSCHNI'^UES 


I. Introduction 

Recently there has been strong research interest on high 
resolution spectral analysis techniques. This is an important 
area of defense research because of the numerous applications to 
radar, sonar, and geophysical areas of defense interest. An 
excellent publication is the Proceedings of the 1976 and 1979 RA^C 
Spectrum Estimation ii/orkshop. Maximum entropy spectral analysis 
is one of a number of high resolution spectral analysis techniques. 
The impact of the Burg's maximum entropy spectral analysis method 
is far more significant than the technique itself. Thus in this 
report we present not only the bibliography of the maximum entropy 
methods in one and two spatial dimensions but also a number of 
related methods of high resolution spectral analysis. One common 
assumption with all these methods is that the data record is short 
and thus the conventional fast Fouriter transform method of spec¬ 
tral analysis is not suitable. Probably because of the snort length 
record, the maximum entropy spectral computation is fairly sensitive 
c the presence of noise. In the following sections, references are 
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